Electron interactions in graphene in a strong magnetic field 
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Graphene in the quantum Hall regime exhibits a multi-component structure due to the electronic 
spin and chirality degrees of freedom. While the applied field breaks the spin symmetry explicitly, 
we show that the fate of the chirality SU(2) symmetry is more involved: the leading symmetry- 
breaking terms differ in origin when the Hamiltonian is projected onto the central (n = 0) rather 
than any other Landau levels. Our description at the lattice level leads to a Harper equation; in 
its continuum limit, the ratio of lattice constant a and magnetic length Ib assumes the role of a 
small control parameter in different guises. The leading symmetry-breaking terms are lattice effects, 
algebraically small in a/Is- We analyze the Haldane pseudopotentials for graphene, and evaluate 
the easy-plane anisotropy of the graphene ferromagnet. 

PACS numbers: 73.43.-f, 71.10.-w, 81.05.Uw 



Introduction: The recent discovery 1,2 of an integer 
quantum Hall (QH) effect in a two-dimensional (2D) 
sheet of graphite, known as graphene, has triggered 
an avalanche of activity, including on the theory side 
studies of the transport properties of relativistic Dirac 
particles* 3 *^ the analysis of edge states&££ and shot 
noised as well as Berry phases in bilayers. 10 

In a simple model, electrons in graphene can be treated 
as hopping on a honeycomb lattice^i^ Perhaps the most 
salient feature of this problem is the existence in the band 
structure of a pair of Dirac points with a linear ('relativis- 
tic') energy-momentum relationship. These points are lo- 
cated at the two inequivalent corners of the Brillouin zone 
(labelled by K and K' = —K), endowing graphene with a 
multi-component structure analogous to the well-studied 
examples of bilayer systems, the multi- valley structure of 
Silicon and of course the simple spin degree of freedom. 
These degrees of freedom can be thought of as SU(2) (or 
higher symmetry) pseudospins. The resulting Hamilto- 
nian typically contains symmetric terms as well as ones 
which lower the symmetry, such as Zeeman (spin) or ca- 
pacitance (bilayer )i£ energies. 

Here, we argue that graphene may be viewed as a fur- 
ther type of multi-component system. Its internal de- 
gree of freedom can be thought of as a chiralityii^ the 
wavevectors K and K' encode the (anti)clockwise vari- 
ation of the phase of the electronic wavefunction on the 
three sites neighboring any given site on one sublattice. 

Moreover, we show that for n^O the chiral SU(2) sym- 
metry is reduced to U(l) in graphene, due to backscat- 
tering terms with momentum transfer 2K ~ K which 
provide a coupling of the chirality to the orbital part of 
the wavefunction. On the contrary, in n = 0, the bro- 
ken symmetry may be due to electrostatic ("Hartree") 
effects. Although different in origin, both effects are of 
order 0{a/ls). The distance between neighboring carbon 
atoms a ~ 0.14nm provides an additional lengthscale be- 
sides the magnetic length Ib = y/h/eB — 26nm/v/B[T], 



so that a -C Ib- It is somewhat analogous to the layer 
separation d for bilayers. In graphene for n^O, however, 
only the exchange part of the symmetry-breaking inter- 
action is non-zero, whereas in bilayers the direct term 
encodes the capacitance energy, which can be important 
already for typical values of d/ls ~ 1. 

In the following, we flesh out this picture with a micro- 
scopic calculation starting at the lattice level, in which we 
derive and discuss the effective model for interacting elec- 
trons restricted to a single relativistic LL and compare 
it to the non-relativistic case of electrons in conventional 
semiconductor heterostructures; the difference between 
the two is most significant for n = 1 . The backscattering 
terms are discussed in the case of the QH ferromagnet at 
the filling factors v = 1 of the partially filled LL (for an 
arbitrary LL, we have v = An + D). 

The model: The electron field in graphene may be 
written as a two-spinor whose components, tpa-( r ) = 
exp(icrK • r)x CT (r), are a product of a slowly varying 
part Xo-(r) and a rapidly oscillating plane wave with 
crK = cr(47r/3v / 3a)e I for the Brillouin zone corners K 
and K' (chiralities a = ±1, respectively). The compo- 
nents of each two-spinor field Xo-(r) correspond to the 
two triangular sublattices (labelled by a — ±1) of the 
bipartite honeycomb lattice. In a magnetic field with the 
Landau gauge A = (0, Bx), q v is a good quantum num- 
ber, and we may expand Xa( T ) — ex p(*9i/y)<?a(y)- The 
electron dynamics (in a tight-binding model with nearest- 
neighbor hopping t = 1) is governed by the Harper equa- 
tion 
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where the distances are measured in units of a. In 
order to derive a continuum limit in the presence of 



an unbounded vector potential A, one expands the co- 
sine in Eq. Q in the vicinity of x^ defined as [q y + 
ax ll (q y )/l B }V3/2 = 2irp, where fi is an integer which ef- 
fectively acts as an additional quantum number besides 
the quasimomcntum q y in the first Brillouin zone. 
The continuum limit of Eq. JQ| thus reads 
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where x is now a small deviation from x^. This result 
coincides with the ones obtained by introducing the min- 
imal coupling p — ► p + eA after deriving the B = 
continuum theoryii^i^ Note that the typical extension of 
the wavefunctions along the x-axis is Rl oc y/nlg in the 
n-th LL, and the periodicity of x^ is ~ 1%/a. The over- 
lap between wavefunctions with differing fi is therefore 
exponentially suppressed provided y/n <C Ib/cl- Finally, 
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where the index + represents the K point. In the ex- 
pression for X- ( a t K'), the components of the spinor 
are reversed. Here sgn(n) = {1,0, —1} for n{>, =, <}0, 
respectively. The quantum number n is the index of the 
relativistic LL, and m is associated to the guiding center 
operator, which commutes with the one-particle Hamil- 
tonian. The \n, m) are the usual (non-relativistic) one- 
particle states for a charged particle in a perpendicular 
magnetic field. The c„ !m;0 - are fermionic destruction op- 
erators. 

Projection onto a LL (n ^ 0) of the sublattice densities 
Pa{r) = E<7, CT ' $L(r)Vv,a(r) gives 



where the projected density operators read p a(J (q) = 
Em,m'( m l exp{-i[q+ (a - a')K] ■ R}|m')ct )m>0 .c„, m ' )0 ,. 
The operator R represents the usual guiding center posi- 
tion, and r\ is the cyclotron coordinate, r = R + rj. The 
chirality-dependent form factors F° a (q) read, in terms 
of associated Laguerre polynomials L"(x): 
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The Hamiltonian of interacting electrons in graphene, 
projected onto a single relativistic LL thus reads 



1 
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where the sum over the wave vectors is restricted to the 
first Brillouin zone. Indeed, the potential consists of a 
sum over reciprocal lattice vectors, as the local densities 
[Eq. ©], valid for l B 3> a, are restricted to the hexagonal 
lattice. The chirality-dependent effective interaction 
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is not SU(2) symmetric; however, the symmetry-breaking 
terms are suppressed parametrically in a/Is- To see this, 
consider the different form-factor combinations in the ef- 
fective interaction potential J3J. 

- Terms of the form F£ a (Tq)^' ~°' (±q) and "umk- 
lapp scattering" terms [F°~ a {— q)F£ ,_<T (q)] are expo- 
nentially small in a/Is ■ 

- "Backscattering" [F^ a (-q)F~ (T << T (q)]: one obtains 
v+-+(q) ~ exp(-|q| 2 /2)/|q'±K| ~ exp(-|q| 2 /2)/|K|, 



which is only algebraically small, 
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a/Is, and thus constitutes the leading perturbation to 
the remaining [SU(2) invariant] terms. 

These leading-order terms in the effective interaction 
yield the SU(2) [or SU(4), if the physical spin is also 
taken into account] symmetric Hamiltonian (for n ^ 0) 
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with the graphene form factor 
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and Pcr(q) = /5 CTCr (q). The graphene form factor (JJJ has 
already been written down by Nomura and MacDonald 
in their study of the QH ferromagnetism at v = l. 16 
The leading-order symmetry-breaking correction due to 
backscattering is (with <' _<r (q) = < ,_ °'' _0 '' <r (q)) 
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The central level (n = 0) behaves remarkably differ- 
ently. In this case, the electron chirality a is equivalent 
to the sublattice index a (Eq. [5J , and therefore 
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where q — q x + iq y and K = K x + iK y are written in 
complex notation, and the wave vectors are given in units 
of 1/Ib- F+ + (q) = F n (q), and F~ + (q) is obtained by 
replacing q — > — q in F^ (q). 



with the same form factor J 7 n= o(q) = e~ q ' 4 as for non- 
relativistic electrons in the lowest LL. From an electro- 
static point of view, it may be energetically favorable 
to distribute the electronic density with equal weight on 
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FIG. 1: Effective interaction potentials in real space: compari- 
son between the relativistic (black lines) and non-relativistic (gray 
lines) LLs n = 0,1, and 5. The dashed line represents the inter- 
action potential in n = 0, which is the same for relativistic and 
non-relativistic electrons. 



both sublattices. For n/0, this follows directly from 
Eq. J2J, but in n = 0, an equal- weight superposition 
of a = ±1 is required to distribute the charges homoge- 
neously on both sublattices. Such an electrostatic effect, 
compared to the SU(2) invariant terms, is of the same 
order 0(cl/Ib) as the backscattering term in u ^ 0. 
Effective interaction: By absorbing the form factor 
into the interaction, we define 
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Fig-fflshows the effective interaction potentials © trans- 
formed to real space, for n — 0,1, and 5. At large dis- 
tances, the usual 1/r Coulomb potential is obtained. In- 
terestingly, the shape of the interaction potential for the 
relativistic n = 1 LL in graphcnc is more similar to the 
non-relativistic n = level than to the corresponding 
one n = 1, as may also be seen in a pseudopotential 
expansion, 17 V£ = {2n)- x £ q v„(<z)M9 2 ) exp(-q 2 /2). 
Indeed the ratios V2m+i/V2m+3 for the odd integer pseu- 
dopotentials, which are relevant for the case of polar- 
ized electrons, decrease monotonically both in n — 1 and 
n = for relativistic LLs. These ratios, and the differ- 
ences Vbm-i — V2m+i are bigger in the former case, so 
that - among the polarized states - fractional QH states 
will therefore be more stable in n = 1 than in n — (at 
constant magnetic field). By contrast, candidate chiral- 
ity unpolarized states (such as at v = 2/3) fare better, 
for two reasons: firstly, the fact that the relativistic ef- 
fective potential is more short-ranged in n = than in 
n = 1 leads to Vb (and Vo/Vi) being smaller for n = 1. 
Secondly, the pair of internal SU(2) degrees of freedom 
allow for a smaller unpolarized 'composite Fermi sphere'. 
Numerical results^ show a first-order phase transi- 
tion at V1/V3 ~ 1.3 between the Pfaffian stated at 
v = 1/2 and a charge-density wave^Si and a crossover 
to a composite- fermion Fermi sea when V1/V3 is fur- 
ther increased. The Pfaffian state is absent in n = 0, 
where V® /V$ = 1.6, probably due to an inaccessibly 
small gapJ^ In the relativistic n = 1 LL, one finds 



an even larger ratio T^ 1 1 /V 3 1 = 1.67 so that a Pfaffian 
state is also unlikely to be observed there. Even though 
the ratio V\ jV% — 1.16 in the relativistic n = 2 LL 
is larger than in the corresponding non-relativistic level 
(V1/V3 — 1.14), it is well below the critical ratio, and 
one would thus expect a stripe phase at v = 1/2. 

It is straightforward to check that the difference to the 
non-relativistic case vanishes, in the large-n limit (see 
Fig- n for n = 5), i.e. far from the Dirac points. Re- 
placing L n (q 2 /2)ex.p(—q 2 /A) ~ J (q^2n + 1), the enve- 
lope otJ*(q) ~ [j (qV^r^T) + MqV^TTT)] 2 /i ~ 
jQ(q\/2n) agrees to leading order in n with the non- 
relativistic case. 

QH ferromagnet at v = 1: In recent transport mea- 
surements on a single graphene sheet additional integer 
QH plateaux beyond those corresponding to v = An 
have been observed^ These appear as the first sig- 
nature of electron-electron interactions, and the anal- 
ogy with the non-relativistic case in semiconductor het- 
erostructures hints at a chirality QH ferromagnet. The 
stability of such a state, in the presence of impurities, 
has been investigated by Nomura and MacDonald. 16 We 
now analyze the impact of the backscattering term (JSJ 
on such a ferromagnet for n 7^ 0, within the Hartree- 
Fock (HF) approximation. Following Ref. I13L we con- 



sider the HF trial state \#) = n m ( M r, 
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where we may parametrize u m — cos{9 m /2)e~ l ^ m ' 2 and 
v m = sin((9 m /2)e^ m / 2 , in terms of the real angle fields 
9 m and <J) m , which can be thought of as polar coordinates 
of a vector field n(m). In the case of a SU(2)-symmetric 
repulsive interaction, it has been shown that the trial 
state |\P) minimizes the energy for constant 9 m and <p m , 
thus yielding a simple ferromagnet »ii The backscattering 
term, averaged over this state, is, apart from an unim- 
portant constant C, 

(H bs ) = - Y^ { V ™,rn>,m,m> [n x (m)n x (m') + x -» y] 
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The factor of \q\ 2 in this sum is due to the fact that the 
wavefunctions on the same sublattice, but for different 
chiralities, are orthogonal. A gradient expansion!^ yields 
to lowest order an easy-plane anisotropy A z : 
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This is reminiscent of the bilayer case, where a finite layer 
separation also induces easy-plane ferromagnetism. The 
key differences are: (i) the parameter a/Is ~ 10~ 2 , which 



mimicks the "layer separation" , is tiny for currently ex- 
perimentally accessible magnetic fields. This implies a 
Curie temperature 9 ~ e 2 /ksdB, whereas the crossover 
to easy-plane behavior does not become visible until a 
logarithmically (in a/Is) small energy. As chirality fer- 
romagnetism involves neither electric nor magnetic dipole 
ordering, inter-plane coupling in a multi-layer system will 
be suppressed. This opens the perspective of probing the 
2D behavior for instance in specific-heat measurements, 
(ii) Contrary to the bilayer case and the relativistic n = 
LL, the gap is not due to a charging energy when only one 
layer is filled - there is no contribution to Eq. (|12[1 from 
the direct interaction because v^(q = 0) = [Eq. (JSJ)]. 
(iii) A z is a lattice effect - it vanishes linearly in a as 
the lattice constant tends to zero at fixed Ib ■ It does not 
depend on n, whereas the SU(2) symmetric terms scale 
as e 2 /ey/n in the large-n limit. Note, however, that the 
continuum limit based on the Dirac equation ceases to 
be valid when R^ ~ \/nlg ~ 's/a- 

Comparison with experiment: Zhang et aJ. have 
observed additional Hall plateaux corresponding to v — 
0, ±1 (n = 0) and v = ±4 (n — 1).~ The former pair cor- 
responds to a complete resolution of the fourfold degen- 
eracy of LLs corresponding to different internal (spin and 
chirality) degrees of freedom in n = 0. An explanation of 
this has to consider the size of the disorder broadening 
of the LLs, L, compared to their splitting due to the cost 
of exciting quasiparticles away from integer filling^ An 
experimental estimate yields L ~ f .7meViSi 

Using our above results, we find that these quasipar- 
ticles are Skyrmions for n = 0,1, whose energy cost 
E s k = ^np s is obtained within the non-linear sigma 
modelfM with the help of the stiffness 
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One obtains for the experimentally relevant parameters 



(at 17T with dielectric constant e ~ 5) E s k — 47rp s = 
l l ^F/2e 2 /el B ~ I.8meV (n = I) and, for n = 0, E sk = 
\^r/2e 2 /d B ~ 4meV, both for SU(2) or SU(4)^i 
In addition, there is a contribution from anisotropies, 
mainly the Zeeman effect (about Ez = 0.Ii?[T]meV); the 
chirality-symmetry breaking due to lattice effects, being 
of order ~ 0.05meV, play only a minor role here. 

The activation gap at v = ±4 scales linearly with B, 
indicating a relevant Zeeman effect, and the plateau is 
visible from ~ 17T onwards^ Given the Skyrmions in 
n = are more costly than the sum of E s k and Ez 
inn = 1, this explains why the chirality Landau levels 
are resolved at 17T in n — 0, even without the help of 
an anisotropy field, whereas they remain absent at v — 
±3, ±5 in fields up to 45T. In fact, for n = 1, E s k does 
not reach 4meV for fields below 80T; also, the plateau at 
v = disappears below f fT, where E sk ~ 3meV. 

To summarize, we have analyzed a microscopic model 
for interaction effects in graphene in the QH regime. We 
find corrections to the SU(2) chirality-symmetric model 
to be numerically much smaller than the Zeeman energy 
breaking the SU(2) spin-symmetry. In addition, the effec- 
tive interaction potential differs from the non-relativistic 
case most strongly for small but nonzero n, in particu- 
lar n = 1, which will therefore a good place to look for 
interaction physics different from the GaAs heterostruc- 
ture. Finally, recent experiments suggest the presence of 
chirality ferromagnetism and Skyrmions in graphene. 
Note added: After submission of this manuscript, arti- 
cles of related work appeared, by Alicea and Fisheo^i on 
ferromagnetism at the integer QHE, and by Apalkov and 
Chakraborty2£ on exact diagonalisations in the fractional 
QH regime using the above-mentioned pseudopotentials. 
Acknowledgements: We thank N. Cooper, J.-N. 
Fuchs, D. Huse, P. Lederer, and S. Sondhi for fruitful 
discussions. M.O.G. is greatful for a stimulating interac- 
tion in the LPS' "Journal Club Mesoscopic Physics". 
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